Nonequilibrium dislocation dynamics and instability of driven vortex lattices in two 

dimensions 
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We consider dislocations in a vortex lattice that is driven in a two-dimensional superconductor 
with random impurities. The structure and dynamics of dislocations is studied in this genuine 
nonequilibrium situation on the basis of a coarse-grained equation of motion for the displacement 
field. The presence of dislocations leads to a characteristic anisotropic distortion of the vortex density 
that is controlled by a Kardar-Parisi-Zhang nonlinearity in the coarse-grained equation of motion. 
This nonlinearity also implies a screening of the interaction between dislocations and thereby an 
instability of the vortex lattice to the proliferation of free dislocations. 

PACS numbers: 74.60. Ge, 61.72.Bb 



I. INTRODUCTION 



Transport of periodic media such as vortex lattices 
in superconductors and charge-density waves through 
random environments plays a paradigmatic role in con- 
densed matter physics. While the pinning dominated 
low-drive regime exhibiting glassy features has long been 
a subject of extensive research, the nontrivial proper- 
ties of the high-velocity regime were recognized only re- 
cently. The prediction ai a disorder induced nonequi- 
librium phase transition^ from plastic to coherent mo- 
tion of the vortex lattice upon increasing drive, triggered 
extensive studies of rapidly driven disordered lattijCi 
and have atifcrsijCted much recenL-jnterest in theory,El 
experimentsjZrEl and simulations.ll3~lla 

The observed dynamic transition can be described 
qualitatively as dynamic "melting" in analogy to the 
equilibrium melting transition, wiere disorder induces an 
effective "shaking" temperatureld However, this analogy 
cannot be extended to the important questions regard- 
ing the role of disorder in structural transformations of 
periodic structures. The intrinsic nonequilibrium nature 
of the driven state renders the analysis of the structural 
transition specifically challenging. Even the equilibrium 
melting is described theoretically in the two-dimensional 
case only. Within the Kosterlitz-Thouless (KT) theoryEEl 
melting is mediated by the unbinding of dislocation pairs. 
The underlying melting mechanisms are by far less un- 
derstood in higher dimensions. Yet the issue of the sta- 
bility of the ordered phase with respect to the formation 
of topological defects was identified as a key issue for 
the structural transitions, and substantial progress was 
achieved in understanding— the role and contribution of 
disorder in the static case.EZi Whereas in equilibrium the 
criteria for the stability of topological order follow from 
comparing relevant energy scales, the analysis of defect 
nucleation in the nonequilibrium situation is more sub- 
tle since it can no longer rely on energy balance con- 
siderations. In this paper we undertake a study of the 



stability of the topological order of the dynamic state fo- 
cusing specifically on the two-dimensional vortex lattice 
driven through a disordered environment. We describe 
an intrinsic nonequilibrium mechanism giving rise to the 
proliferation of topological defects and therefore the in- 
stability of the driven vortex lattice. 

The outline of this paper is as follows: in Sec. |l| we 
specify the coarse-grained equation of motion for the dis- 
placement field. In Sec. Ill we study the anisotropic 
and subdiffusive nature of the dynamics of a single dis- 
location neglecting the effect of the Kardar-Parisi-Zhang 
(KPZ) nonlinearity, which is discussed further in Sec. 



[V. There we show that the KPZ term leads to a "spi- 
ral" structure of dislocations, screening their long-range 
interaction, and recovering a normal linear diffusion of 
defects. Our results are summarized and discussed in 
Sec. 0. Some technical aspects are deferred to the Ap- 
pendix. 



II. MODEL 

We examine a dilute system of "test" dislocations em- 
bedded into the elastic medium of the two-dimensional 
vortex lattice. The vortices are driven along one of the 
principal lattice directions, the x axis. The dynamics can 
be formulated in terms of a Langevin equation of motion 
for the displacement u of vortices from their perfect lat- 
tice positions that move with average velocity v.t3 The 
motion is governed by a competition of elastic interac- 
tions between vortices, thermal noise, and pining forces. 
Pinning is described by a potential V with local correla- 
tion V{IV)V{0) — Ao(5(R), where the 6 function is sup- 
posed to have a width of the order of the superconducting 
coherence length ^. 

Since the unbinding of dislocations is controlled by 
their interaction on large scales it is legitimate to use 
a coarse-grained] description of the vortex lattice. It was 
derived recentlyo'CTEl that on large spatial scales the equa- 
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tion of motion assumes the form 

ijii = cV^u + SF + C + x^xu + -^(Vu)^ + /(r + vi). (1) 

For our purposes it is sufficient to retain only the com- 
ponent of the displacement field parallel to the velocity. 
Although the other components also experience fluctua- 
tions on large scales, they can only further increase the 
instability of the lattice, which we find below to occur 
even in the absence of transverse displacements. Equa- 
tion (|^) is written in the frame moving with the average 
velocity of the vortices, where each vortex has a vanishing 
average displacement around its average position r. i] is 
the vortex friction coefficient. For simplicity we consider 
the elastic interaction as uniform, i.e., we use only one 
elastic constant c, ignoring a distinction between shear 
and compression moduli and additional anisotropic cor- 
rections obtained from coarse graining. is a thermal 
noise with temperature T. Although 77, c and T are 
renormalized under coarse graining, the corrections are 
small in comparison to the original values for sufficiently 
large drift velocities. The stress coefficierLL.x and the 
Kardar-Parisi-Zhang (KPZ) nonlinearity Ajl3 which are 
absent in the bare equation of motion, are generated un- 
der coarse graining. The pinning force /, which is simply 
a derivative of the pinning potential in the bare equa- 
tion of motion, acquires a random-force character under 
renormalization. In principle, it depends on the displace- 
ment as f{r + vt + u). However, this dependence can be 
neglected after the coarse-graining has been performed 
and / has acquired random-force character. In the limit 
of large drift velocities these parameters take the follow- 
ing values :□ 
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(2b) 

(2c) 

(2d) 
(2e) 



/ and ^ are assumed to be Gaussian distributed with zero 
mean, a is the vortex lattice spacing. The actual driv- 
ing force F = rjv + SF required to achieve the prescribed 
velocity v is determined from the consistency condition 
that the average displacement must vanish. 

The linear stress term in Eq. (|l|) can actually be elim- 
inated by a transformation of the displacement u{r, t) — 
u{r— {x/'i])ti.,t) to a new frame masang with velocity 
V — V — x/v the laboratory frame.Ej The coordinates 
are related by f = r -|- (tx/r])ic, where x denotes the unit 
vector along the velocity direction. The following analy- 
sis is based on the transformed equation. For simplicity 
of notation we subsequently drop the tilde identifying 
transformed quantities. 



III. LINEAR PROBLEM (A = 0) 

We first examine the model in the absence of the KPZ 
nonlinearity to set the stage for its subsequent inclusion. 
In the above model dislocations are incorporated as a 
discontinuity in the displacement field of amplitude a, 
the vortex spacing. Due to the restriction of our con- 
sideration to the displacement component parallel to the 
velocity, the only possible orientations of Burgers vectors 
are parallel or antiparallel to the velocity. 

In order to derive the dynamic response of a dislocation 
to the fluctuations of the elastic medium it is convenient 
to split the displacement field into a topologically non- 
trivial (multi-valued) part uq and a single-valued part 



u{r, t) = -ito(r - rd(t)) + M^(r, t). 



(3a) 
(3b) 



The dislocation position is rd(t) and (/5(r) is the angle en- 
closed between the x axis and r. The problem at hand is 
then to derive an equation of motion for and from 
Eq. (|). To this end we consider for the moment (t) as 
given and find 

■qur^ - cV^Ur^ = C + /(r + vi) - -qid ■ Vuo(r - i"d)- (4) 

Since this equation is linear in one can easily calculate 
Ur^ for given C,, /, and dislocation trajectory r^. Eventu- 
ally we determine the dislocation velocity in response 
to the forces C and / from a condition of local equilibrium 
for the dislocation core. 



Vu^(rd(i),<) = 0. 



(5) 



This equation is valid for "slow" changes of such that 
the displacement singularity "instantly" moves to a point 
where it is no longer subject to forces because its envi- 
ronment is locally homogeneous, can be directly 
obtained from Eq. (Q) that has to be imposed with a 
short-scale cutoff of order a. 

Supposing that the fluctuations of the elastic medium 
are weak and the dislocation displacement is slow, the 
equation of motion for the dislocation can be written in 
the form 



- iwr7d(w)rd(w) = fd(w) 
with a Peach-Kochler-like force 



(6) 



(7) 



(e is the totally antisymmetric tensor; see the Appendix 
for some intermediate steps). Here u™* is the contribu- 
tion to arising from the external forces C, and / at zero 
dislocation velocity. The dislocation drag coefficient 
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diverges logarithmically for small frequencies because of 
the long-ranged response of the displacement field to 
the dislocation motion. This divergence with vanishing 
frequency is equivalent to a divergence with increasing 
system size at zero frequency. It is characteristic for 
two dimensions and has been found previously for dis- 
locations in_iiprtex lattices retaining both displace.rnent 
componentsEiJ and also in pattern forming systems.E3E3 
The force entering Eq. ^ has zero mean and correla- 
tions 



(9a) 



r ~ ~ 2rjT + gS(uj + k • v) 
Jit, rj'^uj'^ + c^k^ 

where the disorder contribution is proportional io g = 
Aq/^^t^^w^ and ka = X]/3 ^apkp. It is instructive to com- 
pare the relative strength of thermal and pinning contri- 
butions 
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for (10b) 



for \uj\ < — . (10c) 



At low frequency the pinning contribution to the y com- 
ponent is negligible in comparison to the thermal contri- 
bution, ^^y^{uj) ^ ^yy{uj). Howcvcr, the pinning contri- 
bution to the a:-component dominates over the thermal 
contribution, ^^'^{uj) > ^^^^i^j). 

Since fd vanishes on average, the dislocation has mean 
zero velocity in the frame of Eq. (^), i.e., their velocity 
is smaller than that of the vortices by x/v oc v~^. The 
diffusive behavior of dislocations is readily found from 
Eq. (^). In the absence of disorder dislocations actually 
move subdiffusively, 



([rda(i)-rda(0)]2) 



32T 



\t\ 



(11) 



for \t\ 3> rjo? / Attc. Despite of this peculiar dynamic be- 
havior the fluctuation-dissipation theorem holds for dis- 
locations in the absence of disorder, which implies that 
dislocation pairs nevertheless unbind at the KT transi- 
tion temperature. 

The presence of disorder can significantly affect the 
diffusive behavior of the dislocation. For qualitative 
purposes we may still use the linear response approach, 
where the forces acting on the dislocation are evaluated 
at the undisplaced dislocation position. This approach 



^gjj^i^uld provide a good approximation for a considerable 
time interval because of the logarithmic divergence of the 
drag coefficient. From the above equations one finds that 
motion is super diffusive along the x-direction. 



([xd(i)-xd(0)]2 
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1/2 



(12) 



on large time scales \t\ ^ if'v'^T^/cg^, where the pin- 
ning contribution dominates over the thermal contribu- 
tion. The validity of Eq. ( p^ ) is limited by the fact that 
the nonlinear terms in the equation of motion for rd have 
been neglected. This could lead on largest time scales 
to a further renormalization of the dislocation velocity 
V. However, a quantitative calculation of the dislocation 
velocity is beyond the scope of the present paper, since it 
would require the inclusion of additional contributions to 
the equation of motion, such as those due to Peierls bar- 
riers. However, despite of the possibility that Eq. ( |T^ ) 
may not capture the true large-scale behavior one may 
conclude from a comparison of Eq. ( |l^ ) to Eq. ( |ll]) that 
the shaking effects of disorder correspond on large scales 
to an effectively infinite temperature. This result pro- 
vides a first indication for a disorder-driven unbinding of 
dislocations. 



IV. NONLINEAR PROBLEM (A = 0) 

So far the quadratic KPZ nonlinearity has been ex- 
cluded from our analysis and now we address the ques- 
tion how it modifies the above findings. For simplicity, 
we initially drop the random force term and the thermal 
noise to examine the structure and dynamics of single dis- 
locations as well as the interaction between dislocations. 
In this case Eq. (n^) is reduced to the Burgers equation 



SF. 



(13) 



In analogy to the usual pinning problem we consider 
the vortex drift velocity, which enters the last equation 
through A, as prescribed and determine the related force 
contribution SF from the stationarity condition u — 0. 
Topological defects in Eq. (|l3|) were studied recently in 
the context of pattern- forming sy stems, c2l"E3 where they 
constitute the source of "spiral waves." Therefore it is 
possible to carry over part of the previously achieved 
analysis to the present context. 

It is convenient to perform the well-known Hopf-Cole 
transformation that leads to a linear equation for the 
function W = exp{Xu/2c), 



riW ^cV^W - cklW. 



Looking for a stationary solution, the constant 



(14) 



(15) 
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will be determined from the condition that the solution 
is not singular at the core (for the simplest choice ko = 
the solution W exhibits a non- physical singularity at 
r = 0). Note, that the stationarity of u implies that A 
and SF must have opposite sign. 

Equation ( |l3| ) possesses a solution of the form 



u(r) = —f{r) + fi{r), 
Ztt 



(16) 



where fi(r) is a rotation-symmetric contribution induced 
by the KPZ nonlinearity. This solution still repre- 
sents a topological defect with the characteristic multi- 
valuedness of the displacement u{ip + 2Tnn) — u{ip) + am 
(our explicit considerations apply to m = 1 only, the gen- 
eralization to integer m is straightforward) . The angular 
and radial dependences of W factorize, 



Wo{r,ip) — exp[(aA/47rc)(/3] w{r) 



(17) 



and for a stationary solution of Eq. ( p^ the radial con- 
tribution w(r) = exp[A/x(r)/2c] has to satisfy a modified 
Bessel equation 



„2 1 o a ,9 
a^w H — drW H — -w — KqW, 



(18) 



where a = aX/Anc. Thus the solution to Eq. ( [l8|) is a 
modified Bie^el function z/;(r) = Kia{kor) with imagi- 
nary index.E3 This solution has two characteristic length 
scales separating three regions. On large length scales 
r ^ Lg = k^^ this solution decays exponentially, w(r) oc 
exp[— fcor]/r^/^. For r <C kQ^ and a ^ 1 an expansion 
of the Bessel function yields w{r) w sin[a(ln(2/fcof') —7)] 
with Euler's constant 7. Thus w{r) assumes a maximum 
at a second characteristic scale 



ro = fco 



(19) 



For r 4C ro the solution becomes strongly oscillating, 
which is unphysical. It is important to keep in mind that 
the equation of motion for the displacement field is valid 
only on scales larger than a cutoff scale Rc of the order 
of the vortex spacing. Therefore the oscillatory behavior 
of w at small scales is an artifact of Eq. (|l^) that does 
not account for the dislocation core structure. The above 
solution is physically meaningful only as long as fi{r) or 
w{r) depend monotonously on the distance, i.e., the out- 
ermost maximum of the solution found above has to be 
identified with the core radius. Thus Eq. (|l9|) actually 
determines fco: <p*iequivalently the screening length Lg as 
a function of q;,E2I 



Ls = - Rc e'^/^" ~ a e 
3 _ AqC 



(v/vof 



(20a) 
(20b) 



For small a, i.e., large velocity v of the vortex lattice, 
Lg is exponentially large. Since for r — > 00 the solu- 
tion w ^ exp[— fcorj/r^/^ has exponential asymptotics, 



the displacement u « (2c/A) In w « — (2c/A)A:or increases 
proportional to the distance from the dislocation core. 

The exponential decay of the function K readily im- 
plies exponential screening of the interaction between sev- 
eral dislocations at positions with distances large com- 
pared to Ls- For pattern forming systemSj— this screen- 
ing was a subject of intensive investigationsEjcIl Indeed, 
since Eq. (14) is linear in W, the multi-dislocation so- 
lution can be approximated as a linear superposition of 
individual solutions, 



W{r) = J2Wo{\r-ri\) 



(21) 



where the correction li is introduced in order to fix "topo- 
logical conditions" imposed by the field u. For a well- 
separated ensemble of dislocations the overlap between 
the individual contributions to W is exponentially small. 
From the velocity of a dislocation under the influence 
of other dislocations one can actually obtain the dis- 
location interaction. The resulting interaction between 
two dislocations i and j decays as exp[— 2fco|ri — rj|] for 
A:o|i"i ~ Tjl ^ 1 (for fco|ri — rj| <-|L-there is a crossover to 
the usual power-like behavior). cSEil 

The structure of such a "spiral" dislocation is illus- 
trated schematically in Fig. [l. Note that the lattice is 
compressed in front of the dislocation, whereas it is di- 
luted behind the dislocation. This effect is independent of 
the orientation of the Burgers vector parallel or antiparal- 
lel with the velocity. Given the fact that A is positive this 
effect can be easily understood from Eq. ( p^ ) , if one con- 
siders the KPZ nonlinearity as perturbation to the usual 
dislocation: the displacement strains are largest close to 
the dislocation core. Therefore the vortices close to the 
core experience a force that drives them further in the 
direction of the velocity until the elastic forces establish 
force equilibrium in a configuration, where all vortices 
move with the same velocity v. 

Due to the presence of the KPZ term in the equation of 
motion (|l3|) the force required to achieve a vortex motion 
with velocity v is reduced: SF ^ F — rjv — ~{2c/X)k'^ < 
0. However, this contribution of the KPZ term, which is 
generated by disorder, actually represents only a further 
correction to the immediate pinning force ^p), which is 
positive. For large v the correction due to the KPZ term 
is negligible in comparison to the pinning force, since 
the former contribution decays exponentially for large v 
whereas the latter decays only algebraically. Consider- 
ing both contributions consistently together, one finds 
F > rjv, i.e., even in the presence of the KPZ nonlinear- 
ity disorder actually slows down the vortex motion. 

The structure of the vortex lattice in the presence of 
a dislocation/antidislocation pair is illustrated in Fig. ||. 
This structure was calculated numerically directly from 
Eq. ( p^ using a lattice domain with periodic boundary 
conditions. In Fig. ^ two such domains are reproduced. 
The inhomogeneity of the density is apparent and the 
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density can change along narrow boundaries that corre- 
spond to shock wave fronts.E3 

Let us now include the effect of weak thermal noise and 
disorder in the dynamics of a single dislocation. After the 
Hopf-Cole transformation Eq. (|l|) assumes the form 

TjW = cV^W-cklW+—[C + fir + vt)]W. (22) 

2c 

To derive the resulting equation of motion for a disloca- 
tion we apply a perturbation technique]^ where the force 
fields of thermal noise and disorder are projected on the 
translation modes Wx^y of the unperturbed Eq. ( p2|) (see 
also the Appendix). Since the translation mode is simply 
Wx,y = dx,yWo, the equation of motion for the disloca- 
tion position ro is of the form 

xWo(r')VWo(r') (23) 
with an effective drag coefficient ijd of the form 

m = lj d'r'[VWo{r')r. (24) 



Because of the exponential localization of the function 
Wo, the drag coefficient rjd is finite, in contrast to the 
linear case, where it diverges logarithmically with the 
volume. From Eq. ( p^ ) we obtain the mean squared dis- 
placement [abbreviating C*^"^ = C(r^"-* , i^"-*) and /'"^ = 
/(r(") -ro(t("))-hvt("))] 

- /"* dt' f dt" f d\'d\"{[C + f][C" + /"]) 



'ic^Vd Jo JO J 
xWo(r')VWo(r")Wo(r")VWo(r") 



(25) 



Let us first consider only the effect of thermal noise 
Averaging Eq. (|5|) and utilizing the exponential de- 
cay of W at large r, we readily obtain a normal diffusive 
behavior, 



(26) 



Thus, the exponential screening of the dislocation field 
due to the KPZ nonlinearity is responsible for the tran- 
sition from subdiffusion [Eq. (|ll|)] to normal diffusion. 

To include the effect of quenched disorder we have to 
perform an additional disorder averaging of Eq. (25). 
Thereby we have to assume that the position of the dis- 
location is not correlated with the disorder, which is true 
at large enough drift velocity where the dislocation moves 
with the vortices. As in the case of purely thermal noise 
we obtain normal diffusion of the dislocation due to the 
exponential localization of VF, in contrast to superdiffu- 
sion in the linear case. 



V. CONCLUSIONS 



The most important implication of the exponential 
screening of the dislocation interaction is that even ar- 
bitrarily weak thermal noise or random force result in 
an unbinding of dislocations (see Ref. psf ). This means 
that the corresponding Kosterlitz-Thouless temperature 
is zero in this situation, and that the topological or- 
der of the vortex lattice is always destroyed on largest 
scales. While previous arguments in favor of such an in- 
stability were based on scaling argumentsEI and numer- 
ical simulationalj we have presented here an intrinsic 
nonequilibrium mechanism. This instability mechanism 
is primarily based on the presence of the KPZ nonlin- 
earity in the equation of motion, which is generated by 
coarse graining the equation of motion of a driven vortex 
lattice. 

At this point we would like to recall that our analysis 
was based on the simplified model with only one displace- 
ment component. Fluctuations of the second component 
will probably lead to a significant further reduction of 
the stability of the vortex lattice, i.e., the true screening 
length could possibly be smaller than the value for Lg 
calculated above. 

It is instructive to compare the dislocation screening 
length Lg to the crossover scale where the KPZ nonlin- 
earity becomes relevant for thermal fluctuations of the 
displacement field in the absence oL-dislocations. The 
latter scale is = exp(87rc^/77TA^) i.e., it is also ex- 
ponentially large for small A. However, the functional de- 
pendence of A is different and L^ <C for small A. [Note, 
that is defined only for finite temperature while Lg is 
well-defined even for zero temperature.] Therefore, dislo- 
cations become important for the structure of the system 
before it can show a crossover to the strong-coupling be- 
havior of the elastic system. 

In a remote analogy the interaction of dislocation pairs 
in a driven vortex lattice can be compared to the inter- 
action of vortex pairs in a superconducting film. In the 
former case screening is induced by disorder, which in the 
latter case is provided by magnetic fields. In the absence 
of screening both systems would perform a Kosterlitz- 
Thouless transition. Screening suppresses this transition 
as a large-scale phenomenon. However, for large enough 
screening lengths this crossover can still be well pro- 
nounced. The screening length of the vortex interaction 
in films is known to be macroscopic. The screening length 
of dislocations Lg as found above will also be macroscopic 
for V vq because of its exponential increase with the 
vortex drift velocity. It is therefore possible that the 
lattice may be well ordered on experimentally relevant 
scales and that a dynamic melting or freezing can still be 
found as a pronounced crossover. 

In the nonequilibrium mechanism for the dislocation 
unbinding examined above the KPZ nonlinearity played 
a major role. It is worthwhile to point out that even if 
dislocations are artificially suppressed in a purely elastic 
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approach, this nonhnearity can induce dynamic ixansi- 
tions between "rough" and "flat" sHding phases. u The 
above analysis demonstrates that even the "flat" sliding 
phase is actually penetrated by free dislocations on length 
scales beyond Lg- Therefore one might suspect that their 
presence modifies the nature of this roughncning transi- 
tion or even blurs the transition, turning it into a mere 
crossover. Again, this crossover may yet be observable 
because Ls grows exponentially for large drift velocities. 

To conclude, we have found an instability of the sim- 
plified one-component model of a drifting vortex lattice 
to proliferation of free defects. As already mentioned 
above, the presence of a second displacement component 
can only lead to a further increase of fluctuations and 
to a reduction of the instability. We have presented an 
intrinsic nonequilibrium mechanism for such an instabil- 
ity based on the presence of the KPZ nonhnearity in the 
coarse-grained equation of motion for the displacement 
field of a vortex lattice driven in a disordered environ- 
ment. The instability is due to an exponential screening 
of the dislocation interaction on a scale that increases 
exponentially with the drift velocity of the vortex lattice. 
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S. acknowledges support from the Deutsche Forschungs- 
gemeinschaft under Project No. SFB341 and Grant No. 
SCHE/513/2-1. 



APPENDIX: 

In this appendix we give some intermediate steps of 
the calculation leading to equation of motion for the 
dislocation in the linear medium. Furthermore we relate 
this calculation to a projection of the force field onto the 
translation modes of the dislocation, as used in Sec. IV. 

For given rd Eq. can be solved for u^. Equation 
(PI) then implies 



= ea/3V/3U^(rd(i),t) 
= ea,3 J <er'dt'VpG{YA{t) -v',t~ t') 



X [C(r', t') + /(r' + vt') - 7^va{t') ■ S7uo{r' - r^{t'))] 

where G is the Greens function for u^. In linear response, 
^, /, and hence rd are small and one may replace rd(t') 
by rd(i) in the previous Equation. This equation has to 
be understood as force balance between the contribution 
fd arising from the external force fields C and / and the 
friction force r7di"d- The separation of these contribution 
leads to Eq. (|). 

Let us now briefly discuss the connection between the 
approaches used in Sees. Ill and IV to calculate the dis- 
location dynamics. The subdiffusive behavior of disloca- 
tions in the linear system can be easily obtained from the 



equation of motion, taking into account that in the linear 
system the translation mode is Uxy = dx^y^Q ^ 1/r. 
Thus, the effective mass diverges as Inr. For the mean- 
squared displacement we obtain (dropping the disorder 
term) 

lvl{t)) = — f f [ <fir'(fr"dt'dt" 
Vd Jo Jq J 

x(C'C"Vuo(r')V7.o(r")) 

instead Eq. (|25|). After averaging we obtain (r§(t)) ~ 
t/\nr « t/ \nt, which coincides with Eq. {Kw- 
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FIG. 1. Schematic illustration of the structure of a single 
dislocation (marked by the symbol _L) in the vortex lattice 
driven along the x axis (lengths are in units of a). The den- 
sity is increased/reduced in directions in front of/behind the 
dislocation. 




FIG. 2. Schematic illustration of a dislocation pair (re- 
peating two unit cells with periodic boundary conditions) in 
the driven vortex lattice. The vortex density changes quite 
abruptly along pronounced shock wave fronts. The vortex 
lattice is shown here in a square lattice to make the local 
rotations and compression of the unit cell graphically more 
transparent although in reahty the vortex lattice is hexago- 
nal. 
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